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Abstract 

This article combines two talks given by the authors and is based on works 
done in collaboration with G.E. Brown and D.-P. Min on kaon condensation 
in dense baryonic medium treated in chiral perturbation theory using heavy- 
baryon formalism. It contains, in addition to what was recently published, as- 
trophysical backgrounds for kaon condensation discussed by Brown and Bethe, 
a discussion on a renormalization-group analysis to meson condensation worked 
out together with H.K. Lee and S.-J. Sin, and the recent results of K. M. West- 
erberg in the bound-state approach to the Skyrme model. Negatively charged 
kaons are predicted to condense at a critical density 2 <^ p/po Si 4, in the range 
to allow the intriguing new phenomena predicted by Brown and Bethe to take 
place in compact star matter. 



* Based on talks given by CHL and MR at International Workshop on Nuclear and Particle 
Physics, "Chiral Dynamics in Hadrons and Nuclei", Feb. 6 ~ Feb. 10, 1995, Seoul Nat'l University, 
Seoul, Korea 
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1 Motivation 



Recent work by Bethe and Brown [p]] on the maximum mass of stable compact stars - 
called "neutron stars" in the past but more appropriately "nuclear (or nucleon) stars" 
- suggest that the nuclear equation of state (EOS) in the interior of compact stars 
must be considerably softened at densities a few times the nuclear matter density 
Po by one or several hadronic phase transitions. It is now fairly clear that neither 
pion condensation nor quark matter will figure at a density low enough to be relevant 
to the star matter although the issue is not yet completely settled. As Bethe and 
Brown suggest, kaon condensation could however take place at a density 3-4 times the 
normal matter density and hence play an important role in explaining the remarkably 
narrow range of compact star masses observed in nature ||. 

The aim of this talk is to describe a higher-order chiral perturbation calculation 
that predicts the critical density for kaon condensation. The strategy is to take up 
what Kaplan and Nelson|| started, namely chiral perturbation theory (%PT). Kaplan 
and Nelson predicted in tree order of that kaons condense in neutron matter at 
P ^ 3po- Our calculation goes to next-to-next-to-leading (NNL) order. It turns out 
that the calculation confirms the Kaplan-Nelson prediction although in the process 
new and interesting physical elements are uncovered. Our result is that for reasonable 
ranges of parameters involved, the critical density comes out to be 

2 < p/po < 4. (1.1) 

This is the range of densities relevant to the Bethe-Brown scenario for the formation 
of light-mass black holes for stars that exceed the critical mass of M ~ 1.56M . Their 
arguments extend the estimated range of main sequence star masses, for which stars 
go into black holes, down to ~ 18M Q . 

There are two situations where the production of kaons brings out interesting 
physics. One is their properties in relativistic heavy-ion collisions that involve tem- 
perature. Here kaon condensation is not directly relevant but the mechanism that 
triggers kaon condensation in the relevant situation has intriguing consequences on 
the properties of kaons observed in heavy-ion experiments. This is discussed in a 
recent review^] and will not be discussed here. What we are interested in is what 
kaons do in cold dense matter appropriate to compact objects that result from the 
collapse of massive stars. 

In stellar collapse, as matter density p increases, the electron chemical potential 
p e (determined by the chemical potentials of neutrons and protons in the system 
together with charge neutrality) increases, reaching several hundreds of MeV. If the 
electron chemical potential reaches the "effective mass" of a meson $, m$, then the 
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electron can "decay" into a $ as|J 

e~ -> $~ + v e . (1.2) 

In nature, the only low-mass bosons are the pseudo-Goldstone bosons $~ = 7r~, If - . 
While lowest in mass, the pions do not seem to play an important role, so the next 
possible boson is the kaon with its mass ~ 500 MeV in free space. The electron 
chemical potential cannot reach this high, so on-shell kaons cannot be produced by 
this process. However as will be described below, the kaon in medium can undergo 
a mass shift due to density-dependent renormalization. As the \x e increases and the 
effective kaon mass M* K decreases as p increases, the process (|1.2j ) can occur at some 
density p c . Kaons so produced will bose-condense at that density p c . Whether or not 
this will occur then depends on whether or not will decrease enough in density 
so that it meets p e . Such a condensation will be of physical interest if the critical 
density is low enough and the energy gain is high enough. This is the possibility we 
shall address below. 



2 Maximum Neutron Star Mass 

2.1 Stellar Death Function 

That main sequence stars of mass > 25 — 3OM must end up in black holes without 
producing nucleosynthesis, i.e., without returning matter to the galaxy, is required 
by the observed abundances of elements^, [7J. Maeder's argument is based on the 
measurement of AY/AZ, the ratio of helium abundance to that of metals, in low- 
met allicity extragalactic Hjj regions, especially irregular dwarf galaxies. The ratio 
can be measured with good accuracy ||, 

^ = 4 ± 1.3. (2.3) 

If all stable stars of mass up to ~ 100M Q were to explode, returning matter to the 
galaxy, this ratio would lie between 1 and 2. Helium is produced chiefly by relatively 
light stars, metals by heavy stars, so that cutting off the production by the heavy 
stars going directly into black holes without nucleosynthesis increases the AY/AZ. 
Using the standard initial mass function for stars, 

dN/dM = M~ (1+x) (2.4) 

with x = 1.35,1.70, Maederjf] found that Pagel's[|J measurement on AY/AZ was 
best reproduced by a cutoff of nucleosynthesis at a main sequence stellar mass of 
~ 22.5M Q as shown in Fig. [I]. 



3 



AZ 8 



6 



2 







1 1 

- 


J 1 1 1 1 1 1 1 1 1 1 1 

□ : x = 1.70 - 
■ : x = 1.35 - 










*.■.■.■.■.■.■.■.■.-.-.-. Observed *.*.*.*.• 






I I I I I I I I I I I I l_ 

20 40 60 80 100 120 

M BH 



Figure 1: Values of the ratio AY/AZ of the relative helium to metal enrichments for 
different values of x as defined in eq.( [2.4|) as a function of Mbh- The data points cor- 
respond to initial metallicity Z=0.001. The observed range is indicated by shadings; 
the range 4 to 5 is the preferred one.[0] 

There is considerable uncertainty in the initial mass function, as noted by 
Maederf7j, so that this limit could easily be ~ 3OM or even higher. Brown and 
Bethe estimated 3OM as the cutoff for stars to drop directly into black holes 
without nucleosynthesis. 

What about stars with mass 2OM < M < 3OM ? Recent observation on 
SN1987A, whose progenitor mass is ~ 18 ± 2M , gives us an insight. Based on the 
empirical analysis, Brown and Bethe also argued that a large range of stars below 
this mass, down to ~ 18M , can first accomplish nucleosynthesis and then collapse 
into black holes. 

2.2 SN1987A: Neutron Star or Black Hole? 

SN1987A(February 23, 1987) in the Large Magellanic Cloud is the nearest and bright- 
est supernova to be observed since SN1604AD (Kepler), and certainly the most im- 
portant supernova since SN1054AD, the progenitor of Crab Nebula. Because of its 
brightness and proximity, it will be possible to observe SN1987A for many years as 
it expands to reveal its inner secrets. In contrast, typical supernova, of which some 
20-30 are observed each year, are some 1000 times further and 10 5 times fainter, and 
so become lost in their host galaxy within a year or two. Moreover, SN1987A has 
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Figure 2: The energy of the neutrinos detected by the Kamiokande and 1MB detectors 
is plotted as a single neutrino pulse. Most neutrinos arrived within the first second 
or two but there was a lower energy tail to the pulse which lasted more than 10 sec. 



been observed at every wavelength band of the electromagnetic spectrum, from radio 
to 7-rays. 

In supernova theory, the neutron star is followed by bursts of huge amount 
of neutrinos. In the case of SN1987A, neutrinos were detected by the 1MB and 
Kamiokande detector about three hours before the optical burst. The total energy 
emitted in i? e 's is E u « 3 x 10 52 erg / s, while the decay time scale of burst is At ~ 10 sec. 
These are quite consistent with the expected values. Hence, firstly, it is believed that 
the neutron star was formed in SN1987A, even though there exist time gap of about 
7 seconds between the eighth and ninth Kamiokande events, the gap being followed 
by another three events as shown in Fig. |2][|9], [K| . 

From the observation of radioactive decay of 56 Ni and 56 Co, 

56 Ni 28 (r 1/2 = 6 days) — > 56 Co 27 + e + + 7 + u e 
56 Co 27 (r 1/2 = 77days) — > 56 Fe 26 + e+ + 7 + v e , (2.5) 

the mass of Ni ejecta of SN1987A is known to be O.O75M . Combining this result with 
the observed energy of SN1987A, E = 1.4 ± 0.4/oe (where foe stands for 10 51 erg), 
Brown and Bethejl| obtain the range of the core mass of SN1987A 

M = 1.535 ± 0.02M Q . (2.6) 
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However, from the mass of the Hulse-Tayler pulsar, the lower limit of the compact core 
is known to be at least 1.44M . Consequently, the compact core mass of SN1987A 
must range 

1.44M© < M core < 1.56M . (2.7) 

Astronomers have been searching for a pulsar in the center of SN1987A rem- 
nants after the explosion. To see the explicit signal of a pulsar, one must wait until 
the remnants are transparent. If the pulsar is really formed in SN1987A, the X-ray 
signal should be detected within a few years after the explosion. The hypercritical 
accretion could hide the compact object for only ~ 1 year, but after this time, the 
1987A should be observed with a luminosity L = 4 x 10 38 ergs/s, but the present 
light curve is lower by two orders, L ~ 4 x 10 36 ergs/s|l3]]. This compact object, after 



being a proto-neutron star for at least 10 seconds, appears to have collapsed into a 
low-mass black hole. The core of SN1987A may have become a neutron star followed 
by neutrino emission, and later changed into low-mass black hole with mass about 
1.5M . 

A possible scenario was proposed by Bethe. According to his arguments, a vig- 
orous convection is produced in the supernova shock. But after about 2 seconds, the 
convection stops as heat is no longer supplied by neutrinos. As a result, a substantial 
fraction of the previously convecting material falls into the neutron star at the cen- 
ter. Bethe estimated this fraction to be about 10% of the mass in the shock wave, 
or about 0.04M Q . If the neutron- star mass were close to M max , the added 0.04M Q 
could push the neutron star over the limit and make it collapse into a black hole. If 
this happened in SN1987A, it must have been more than 12 seconds after the first 
collapse, since at 12 sec, a neutrino was still observed at Kamiokande II. A softened 
equation of state might be associated with a delayed collapse of the young neutron 
star into a black hole. 

From standard evolutionary analysis, it is believed that the neutron star is 
formed if the core mass lies between the maximum neutron star mass(M ma2 .) and the 
Chandrasekhar mass,0 

M CH = 5.76Y C 2 M & , (2.8) 
where Y e is the electron fraction per baryon. In stellar collapse, Y e is 0.43 ~ 0.50. 



' In the early 1930's, Chandrasekhar set a limit to the size of white dwarf. No carbon-rich white 
dwarf can support its weight if it is greater than about 1.4 times solar mass Mq. A massive star 
with its final mass after exhausting its fuel is greater than Mqh collapses and the collapse turns 
into an explosion. 
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Because of the thermal pressure, the evolutionary lower bound of the neutron star is 
about (1.10 ~ 1.15) x M CH ~ (1.2 ~ 1.4)M . 

According to Brown and Bet he [0], in the binary pulsar evolution, the accretion 
can proceed at the hypercritical rate 

M > l0 4 M Edd , (2.9) 

where the Eddington limit is 

M Edd = 1.5 x lU-*M Q yr-\ (2.10) 

Hence, if the neutron star mass were determined by the evolutionary scenario, massive 
neutron stars with masses exceeding 1.5M should exist. But as shown below, none 
have been found. 

2.3 Observed Neutron Star Mass 

From radio, optical, x-ray and 7-ray surveys of supernova remnants, Helfand and 
Becker came to the conclusion that nearly half of the supernova in the galaxy leave no 
observable remnants. This is understandable because half of the supernova are Type 
I, which leave no neutron stars. However, as discussed in detail by Van den Bergh et 
al. |M| , the observed samples of supernova in our galaxy are at low galactic latitude, 
so their light is strongly absorbed on its way to Earth. They argued that Type II 
supernova outnumber Type I supernova by a factor of several. Thus the conclusions 
of Helfand and Becker can be understood only if nearly half of the Type II supernova 
explosions do not form neutron stars. The possible candidate is, therefore, a low-mass 
black hole. 

In Fig. H], the measured neutron star masses are plotted. Most of all the neutron 
star masses are below 1.5M Q except for Vela X-l and 4U 1700-37. But by the recent 
analysis]]]], the mass of Vela X-l is believed to be below 1.5M . Further, there is 
an argument that 4U 1700-37 is a low-mass black hole[jL|. If this is the case, it is 
striking that all well measured neutron star masses lie below 1.5M . This calls for 
theoretical arguments to lower the maximum neutron star mass. In Table [3], the 
maximum allowed mass of a neutron star is given for various equations of state |i5|. 

2.4 Lowering Maximum Neutron Star Mass 

Since the pioneering work by Bahcall and Wolf , it has been known that a pion 
condensate, if it exists, will enhance strongly the emissivity of the neutrinos from 
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Figure 3: Neutron star masses in units of solar mass (M ). The empty box at the 
lower end of Vela X-l is given by the recent analysis of Van Kerkwijk.p]] 



Equation of state Maximum mass (M ( 



0; 



7r or K condensate < 1.5 
R 1.6 
BJ 1.9 
TNI 2.0 
TI 2.0 
MF 2.7 

Table 1: The maximum mass of neutron star for various equations of state cited in 
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within the core of a neutron star. These authors considered the reaction 



n + n'-^n + e' + i' (2.11) 
and its inverse reaction, in which the pions were treated as real particles existing 



with a certain probability inside the neutron star. Subsequently Maxwell et al|T7| 
carried out a detailed calculation on the neutrino emissivity in the presence of pion 
condensates. The reaction mechanism may be written symbolically as 

n+ "vr" ^n + e~ +u (2.12) 

where 'V represents the pion condensate built in the neutron quasiparticle states. 



The main conclusion of [0 was that even a small amount of pion condensates will 
cause a dramatic enhancement of neutrino emissivity making the equation of state 
softer than what standard calculations would predict. This soft equation of state 
could then lower the maximum neutron star mass. 

However, the much-discussed P-wave pion condensation is now considered to be 
rather unlikely to take place at a low enough density. Based on renormalization-group 
flow equations, Lee et al||18|| showed that the Yukawa coupling term responsible for 



P-wave condensation, %pf ■ T(^ip, becomes irrelevant after radiative corrections and 
cannot induce instability needed for a phase transition. Furthermore, the axial- vector 
coupling constant g& in nuclear medium is effectively quenched, roughly, down to 
one. ||19|| The quenched axial- vector coupling would then push the critical density to 



pc oc — -, (2.13) 



a higher density (> 5po), 

I 

where most of the approximations associated with effective hadronic Lagrangians 
must have broken down. 

An S-wave pion condensation is also unlikely to occur since chiral symmetry 
protects the pion mass (PCAC). Even if S-wave condensation occurred, the effect 
would be negligible. 

The next candidate process is kaon (K~) condensation. (The kaon mass in free 
space is ~ 495MeV). According to Kaplan and Nelson ||, the attraction needed for 
kaon condensation comes mainly from the KN sigma term 



5Mi « -%^iVtjV) + • • • (2.14) 
f 2 



where T,kn comes from the explicit chiral symmetry breaking 

1 



KN - -ym + m s )(N\uu + ss\N). (2.15) 
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Here rh and m s are ~ 5MeV and ~ 150MeV } respectively. The strangeness content 
of the proton is not well-determined, so that gives some uncertainty in the value of 
T, KN . The predicted critical density was found to be in the range 

2.3p < Pc < 3p - (2.16) 



Recently, an improved calculation was made by Thorsson et al.|20|, How- 
ever this calculation was also incomplete since the tree-order Lagrangian they used 
described correctly neither KN scattering nor kaonic atoms. In these talks, we de- 
scribe how this defect is removed in chiral perturbation theory by going to one-loop 
order pi E3, E3, 



3 Baryon Chiral Perturbation Theory 

3.1 Effective Field Theory for Nuclear Matter 

The process we are interested in requires a field theory that can describe simulta- 
neously normal nuclear matter and phase transitions therefrom. The most relevant 
ingredient of QCD that is needed here is spontaneously broken chiral symmetry. We 
are specifically interested in chiral 577(3) x SU(3) symmetry since strangeness is in- 
volved. In order to address the problem, we need to start from a realistic effective 
chiral Lagrangian, obtain a nuclear matter of the right properties from it and then 
determine whether strangeness condensation occurs. 

Unfortunately we do not yet know how to describe nuclear matter starting 
from a chiral Lagrangian. There are various suggestions and one promising one is 
that nuclear matter arises as a solitonic matter from a chiral effective action, a sort 
of chiral liquid []2B| resembling Landau Fermi liquid. The hope is that the resulting 



effective action would look like Walecka's mean-field model. There is as yet no con- 
vincing derivation along this line. However as argued in [|J, there is a compelling 
phenomenological indication that such a Fermi liquid structure can be identified with 
Walecka's mean field model provided that BR scaling is suitably implemented in the 
fluctuations. In the work reported here, we will have to assume that we have a nu- 
clear matter that comes out of an effective chiral action. Given such a ground state 
containing no strange degrees of freedom, we would like to study fluctuations along 
the strangeness direction and determine if instability along that direction develops 
signaling a phase transition. We are therefore assuming that we can get the prop- 
erties of normal nuclear matter from phenomenology, that is, that nuclear matter is 



a Fermi-liquid fixed point [|27j 2q| . In principle, a precise knowledge of this ground 
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state from a chiral effective Lagrangian at a nonperturbative QCD level would allow 
us to determine the coefficients that appear in the effective Lagrangian with which to 
describe fluctuations around the soliton background - i.e., the Fermi liquid -and with 
which we could then compute all nuclear response functions. At present such a deriva- 
tion does not exist. In a recent paper by Brown and one of the authors (BR91)[f29f . 
it is assumed that in medium at a matter density p ~ p , the nuclear effective field 
theory can be written in terms of the medium-dependent coupling constants g* and 
masses of hadrons m* while preserving the free-space structure of a sigma model. This 
leads to the so-called Brown- Rho scaling. In BR91[E3], the nonlinear sigma model 



implemented with trace anomaly of QCD is used to arrive at the scaling law. The 



precise way that this scaling makes sense is elaborated by Adami and Brown |30| and 
in the review (BR94) |3J . There have been numerous papers written with some of the 
essential points of this scaling misinterpreted. 

Given such an effective field theory, we can make a general argument on the sta- 
bility in various flavor directions of nuclear matter at high density. This can be done 



along the line of arguments developed for condensed matter physics by Shankar[27 



and Polchinski |28] using renormalization group flow. We sketch the essential argu- 



ment following Lee, Rho and Sin|TB 



What we are interested in is whether the system in question develops instability 
along the direction of strangeness and if so, by which physical mechanism. This 
analysis will not give us the critical density. The critical density will be calculated by 
using chiral perturbation theory. For this purpose we will focus on the kaon frequency 



near the electron chemical potential. By Baym's theorem |3i]|, one can identify the 
kaon chemical potential associated with charge conservation, /ix, with the electron 
chemical potential, /i e , which we shall simply write p in what follows. This means 
that we will be looking at the vicinity of u ~ fi in the kaon dispersion formula. We 
shall assume that 

(3.17) 

As mentioned, we assume that nucleons in nuclear matter are in Fermi-liquid state 
with the Fermi energy /j,p and the Fermi momentum hp. Define ip as the nucleon field 
fluctuating around the Fermi surface such that the momentum integral has a cut-off 

k F - A N < \k\ < k F + A N . (3.18) 

Kaons can interact with the nucleons through three-point functions of the KNN type 
(Yukawa interaction) and through four-point interactions of the KKNN type. We 
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shall consider S-wave kaon-nucleon interactions, for which the Yukawa interaction can 
be ignored. A generic action involving the nucleon field tp and the kaon field $ can 
then be written, schematically, as 

S = J dwd 3 q$*{u,q)(u - g 2 /W) ®(u,q) - J dud 3 q M^$*$ 
+ J {dujd 3 q) 2 {ded 3 k) 2 h<5>*<5>^^5\u,e,q,k) 

+ J ded 3 k^ (e - e(Jfe))) V + 9 J (ded 3 k) 4 ^^^5 4 (e, k) (3.19) 

where M = (M^ — fi 2 )/2fi and h and g are constants. The four- Fermi interaction with 
the coefficient g stands for Fermi-liquid interactions in nuclear matter. (In nuclear 
matter, one can have four such terms because of the nucleon spin and isospin degrees 
of freedom. We need not specify them for our purpose.) This is a toy action but it is 
generic in that the results of we will obtain below can be put into this form. 

The renormalization group flow of this action can be analyzed in the following 
way. Since we are assuming that nuclear matter is a Fermi-liquid fixed point, fluctu- 
ations in the non-strange direction in the nucleon sector are stable: The four-Fermi 
interaction g is irrelevant or at best marginal. Fluctuations in the strange direction 
involve the kaon field <E>. Suppose we have integrated out all the high-frequency modes 
above the cut-off A measured with respect to \i. We are interested in the stability 
of the system under the renormalization group transformation A — > sA (s < 1) as 
s — > 0. A scaling analysis shows that the interaction term h is irrelevant while the 
"mass term" M is relevant. The renormalization group-flow of the "mass term" and 
the interaction term h can be readily written down and solved jljj (with t — — In s), 

M{t) = (Mo - ^De* + ^e~ at (3.20) 
1 + a 1 + a 

with 

h(t) = h e~ at , h >0 (3.21) 

where D = Hl+ 2 f )a pN > 0, a = A/k F > and a = 1/2. We see from Eq.(CT) 
that as s — > for which h — *■ 0, M changes sign for some (M ,h Q > 0). Thus 
although irrelevant, an attractive interaction ho determines the direction of the mass 
flow whereas it is the "mass term" that drives the system to instability. 

3.2 Chiral Counting 

Armed with the general information on the instability in the strangeness direction, we 
are now able to calculate the critical density in %PT. As mentioned, we are to look at 
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the instability in the kaon direction, so it suffices for us to look at fluctuations around 
the Fermi-liquid state. For this we need an effective chiral Lagrangian involving 
baryons as well as Goldstone bosons. When baryons are present, x?T is not as firmly 



formulated as when they are absent |32|]. The reason is that the baryon mass me is 
~ A x ~ 1 GeV, the chiral symmetry breaking scale. It is more expedient, therefore, 
to redefine the baryon field so as to remove the mass from the baryon propagator 

Bv = jm BT v*xp +B (322) 

where P + = (1 +7 • v)/2 and write the baryon four-momentum 

= ttibv^ + (3.23) 

where is the small residual momentum indicating the baryon being slightly off-shell. 
When acted on by a derivative, the baryon field B v yields a term of O(k). The (octet) 
baryon propagator simplifies in heavy-baryon formalism to i/v ■ k that involves no 
gamma matrices. This simplifies the loop calculation. The spin operator Sfi is defined 
by, v-S v = 0, S 2 V B V = -\B V , {S£, S£} = \(v*v v and [S*, S%] = ie^v a {S v ) p . 

In the baryon rest frame, the spin operator S v reduces to the usual spin operator a/2. 

Chiral perturbation theory in terms of B v and Goldstone bosons (n ■ A/2) is 
known as "heavy-baryon (HB) xPT" [^] . HB^PT consists of making chiral expansion 
in derivatives on Goldstone boson fields, &m/ 'A x , and on baryon fields, Ob /tub, and 
in the quark mass matrix, nAi/A^. In the meson sector, this is just what Gasser and 
Leutwyler did for tttc scattering. In the baryon sector, consistency with this expansion 
requires that the chiral counting be made with B'(- ■ -)B, not with B(- ■ -)B. This 
means that in medium, it is always the baryon density p(r) that comes in and not 
the scalar density p s (r). This point seems to be misunderstood by some workers in 
the field. 

Following Weinberg ||34j|, we organize the chiral expansion in power Q v where Q 
is the characteristic energy /momentum scale we are looking at (Q << A x ) and 

u = A- N n -2C + 2L + Y,^i (3-24) 

i 

with the sum over % running over the vertices that appear in the graph and 

Aj =di + ~rn - 2. (3.25) 

Here u gives the power of small momentum (or energy) for a process involving N n 
nucleon lines, L number of loops, di number of derivatives (or powers of meson mass) 
in the zth vertex, number of nucleon lines entering into ith vertex and C is the 
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number of separate connected pieces of the Feynman graph. Chiral invariance requires 
that Aj > 0, so that the leading power is given by L = 0, v = 4 — Nn — 2C. 

As an example, consider KN scattering. The leading term here is the tree 
graph with v — 1 and with N n — C — 1. The next order terms are v = 2 tree graphs 
with A = 1 that involves two derivatives or one factor of the mass matrix A4. From 
v = 3 on, we have loop graphs contributing together with appropriate counter terms. 

In considering kaon-nuclear interactions as in the case of kaon condensation, 
we need to consider the case with N n > 2 and C > 2. In dealing with many-body 
system, one can simply fix 4 — N n and consider C explicitly. For instance if one 
has two nucleons (for reasons mentioned below, this is sufficient, with multinucleon 
interactions being suppressed), then we have 4 — N n = 2 but C can be 2 or 1, the 
former describing a kaon scattering on a single nucleon with a spectator nucleon 
propagating without interactions and the latter a kaon scattering irreducibly on a 
two-nucleon complex. Thus intrinsic n-nucleon processes are suppressed compared 
with (n — l)-nucleon processes by at least 0(Q 2 ). This observation will be used later 
for arguing that four-Fermi interactions are negligible in kaon condensation. This 



is somewhat like the suppression of three-body nuclear forces |34| and of three-body 



exchange currents ||35|| in chiral Lagrangians. 



3.3 Kaon- Nucleon Scattering 

Given a chiral Lagrangian, we need to first determine the parameters of the La- 
grangian from available phenomenology. This is inevitable in effective field theories. 
We shall first look at kaon-nucleon scattering at low energies. This was done by Lee 
et al. ^4] which we summarize here. We shall compute the scattering amplitude 
to one-loop order and this entails a Lagrangian written to 0(Q 3 ) as one can see from 
the Weinberg counting rule. Instead of writing it out in its full glory, we write it in 
a schematic form as 

C = Y i C i [B v ,U,M\ (3.26) 

i 

where the subscript % stands for v relevant to the KN channel. Here B v stands 
for both octet and decuplet baryons and U the Sugawara form for octet Goldstone 
bosons. For KN scattering in free-space, the Lagrangian is bilinear in the baryon 
field. Details are given in Lee et al. [22, 24[[. Let us specify a few terms in (|3.26|) so as 



to streamline our discussion. Focusing on S-wave scattering, C\ contains the leading 
order term that may be described by the exchange of an uj between kaon and nucleon, 
attractive for K~N and repulsive for K + N and an isovector term corresponding to 
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the exchange of a p meson. These terms are proportional to the kaon frequency u. 
To next order, £2 contains the " KN sigma term" proportional to T>kn/ f 2 where / 
is the pion decay constant and a term proportional to uj 2 which may be saturated 
by decuplet intermediate states. The v = 3 pieces are counter terms that contain 
terms that remove divergences in the loop calculations and finite terms that are to be 
determined from experiments. The complete S-wave scattering amplitudes calculated 
to the NNL order come out to be 
m B 



K±p 

a 



47r/ 2 (m B + M K ) 



TM K + (d s + d v )M 2 K + {(L S + L v ) ± (g s + g v )} M K 



+da Ai , F 



m B 



4nP(m B + M K ) 



T-M K + (d s - d v )M 2 K + {(L s - L v ) ± (g s 



9v)}M- 



(3.27) 



where Mk is the kaon mass, tub the baryon (nucleon) mass, d s is the t-channel 
isoscalar contribution [(9(<5 2 )], and d v is the t-channel isovector one [(9(Q 2 )], both 
coming from £2, L S (L V ) is the finite crossing-even t-channel isoscalar (isovector) finite 
one-loop contribution [0(Q 3 )] having the numerical values 



L S M K w -0.109 fm, L V M K w +0.021 fm 



(3.28) 



and the quantity g s (g v ) is the crossing-odd t-channel isoscalar (isovector) contribution 
[0(Q 3 )] from one-loop plus counter terms in £3. 

To understand the role of the A*, we observe that the measured scattering 
lengths are repulsive in all channels except K~n ||36||. 



K+p 



K+r 



-0.31fm, 
-0.20fm, 



K~p 

a K ~ n 



-0.67 + i0.63fm 
+0.37 + i0.57fm. 



(3.29) 



Although the experimental K~ N scattering lengths are given with error bars, the 
available K + N data are not very well determined. Since both are used in fitting the 
parameters of the Lagrangian, we do not quote the error bars here and shall not use 
them for fine-tuning. For our purpose, we do not need great precision in the data as 
the results are extremely robust against changes in the parameters. The repulsion 
in K~p scattering cannot be explained from Eq.( 3.27| ) without the A* contribution. 
In fact it is well known that the contribution of the A(1405) bound state gives the 



repulsion required to fit empirical data for S-wave K p scattering |22], ^3|, ^4], |3~7fl . As 
mentioned, we may introduce the A* as an elementary field. It takes the form 

9lM 2 K 



6a%* p 



m B 



4nf 2 (m B + M 



k 



m B =f M K - m A * 



(3.30) 
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which is completely determined given experimental data on the coupling g\* and the 
complex mass m\*. 

There are four unknowns d S}V , g Sj „ in ( |3.27| ) which can be determined from four 
experimental (real part of) scattering lengths Eq. ( |3.29| ). The results are 

d s « 0.201fm, d v « 0.013fm, 
g s M K « 0.008fm, g v M K « 0.002fm. (3.31) 

So far, no prediction is made. However given the parameters so fixed, one can 
then go ahead and calculate the S-wave amplitude that enters in kaon condensation. 
This amounts to going off-shell in the u variable, that is, in the kinematics where 
D 7^ Mk- In doing this, one encounters an ambiguity due to the ui dependence of 
the coefficients d which consist of the U KN sigma term" and "u; 2 term" which get 
compounded on-shell into one term. In the calculation reported in Lee et al. [E2L 



we chose to fix the "ur term" by resonance saturation and leave the "sigma term" to be 



fixed by the on-shell data. The predicted off-shell amplitudes |22|, P4j] agree reasonably 
with phenomenologically constructed off-shell amplitudes. All the constants of the 
chiral Lagrangian bilinear in the baryon field are thereby determined to 0(Q 3 ). 



3.4 Four-Fermi Interactions 

In medium, the chiral Lagrangian can have multi-Fermi interactions as a result of 
"mode elimination." Here we consider four-Fermi interactions, ignoring higher-body 
interactions. We shall see that this is justified. 

As stated above, we need to focus on four-Fermi interactions that involve 
strangeness degrees of freedom. Nonstrange four-Fermi interactions are subsumed 
in the Fermi-liquid structure of normal nuclear matter. For S-wave kaon-nuclear 
interactions, we only have the A(1405) to account for. There are only two terms, 

U-fermion = Cf^XTr + ^A^Il (3.32) 

where C*fi T are the dimension —2 (M -2 ) parameters to be fixed empirically and a k 
acts on baryon spinor. We shall now describe how to fix these two parameters from 
kaonic atom data. 

In order to confront kaonic atom data, we need to calculate the kaon self-energy 
n in nuclei. The off-shell amplitude determined above gives the so-called "impulse" 
term 

= -{p P T f l7^)+PnT« e - e n (uj)) (3.33) 
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where T KN is the off-shell S-wave KN transition matrix. (The amplitude T KN 
taken on-shell, i.e., us = Mk, and the scattering length a KN are related by a KN = 
kir(i+M K /m B y ^ KN ■) Medium corrections to this "impulse" term, obtained from one- 
loop graphs by replacing the free-space nucleon propagator by the in-medium propa- 
gator, shall be denoted as 

Pp 8T«-v{u) + p„<^~») . (3.34) 



These two terms ( p.33|) and (|3.34|) are completely determined by the parameters 



fixed above. The new parameters of the four-Fermi interaction come into play in the 
first two self-energy graphs of Fig|4| (the last two graphs do not involve four-Fermi 
interactions but enter at the same order; they are free of unknown parameters), 



n A *M = -^7 



9a* ( w 



f 2 \<jj + niB — m\* 



[C\*Pp (pn + ^P P ) - ^^Pp} 



-gl (——?■ ) (S^(w) + (3.35) 



where g A * is the renormalized KNA* coupling constant determined in Lee et al. f22|, p4"| 
and S^(cj) is a known integral that depends on proton and neutron densities and Mr-. 
Note that while the second term of ( |3.35| ) gives repulsion corresponding to a Pauli 
quenching, the first term can give either attraction or repulsion depending on the sign 
of (C^*[p n + \p P ] — |Ca* Pp)- F° r symmetric nuclear matter, only the combination 
(Cf* — Ca*) enters in the self-energy. This is an important element for kaonic atom. 
The complete self-energy to in-medium two-loop order is then 



n*M = - (p p t£7V) + p n r*; e >)) - { Pp 5r p K N ^{u) + P jr p \ 

+U a ,(oj). (3.36) 



We now turn to fixing the constants of the four-Fermi interactions based on the 
recent analysis of kaonic atoms by Friedman, Gal and Batty [p^j . For later purpose 
we shall parameterize the proton and neutron densities by the proton fraction x and 
the nucleon density u = p/po as 

Pp = xp , Pn = (1 - x)p , P = u Po . (3.37) 



Now Friedman et aZ.[3~8] found from their analysis that the optical potential for the 
K~ in medium has an attraction of the order of 

AV = M K — Mk ~ -(200 ± 20) MeV at u = 0.97 (3.38) 
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Figure 4: The in-medium two-loop kaon self-energy involving A(1405). Figures a and 
b contain the constants of the four-Fermi interaction and figures c and d are Pauli 
corrections 



with 



This implies approximately for x = 1/2 

{cl-cl)f 



20. 



(3.39) 



(3.40) 



Friedman et al. |38| note that their "nominal" optical potential gives an attraction 
of order of 800 MeV when extrapolated to three times the normal density. We show 
in Table |2] what our theory predicts at higher densities than normal.0 At u — 3, the 
net attraction is only about 1.7 times the one at u — 1. 

Equation (|3.35| ) shows that for symmetric nuclear matter (x = 1/2), the combi- 
nation ((7f* + Cj*) does not enter into the self-energy formula. In order to extract it 
as needed for non-symmetric system as in compact star matter, we need information 

* The numerical values in Tables || and || are slightly modified from the previous results in ]25[ 
which had numerically small errors. Specifically, the modifications in STJ^ N , Eq. (G.3) of |24|1 , 



-n N 



M/£f(0) — —z 
2tt z 



d\k\ 



\k\< 



are responsible for the slight changes in the numerical values. We have verified that our new results 
(See Eq. (G.l) of [24]) satisfy the chiral symmetry constraint of MeiBner et al.pjl. 
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u 




Ay 


0.5 


354.9 


-140.1 


1.0 


294.4 


-200.6 


1.5 


249.2 


-245.8 


2.0 


211.7 


-283.3 


2.5 


180.5 


-314.5 


3.0 


153.9 


-341.1 


3.5 


130.8 


-364.2 


4.0 


113.2 


-381.8 



Table 2: K effective mass(M^-) and the attraction (AV = M^ — Mk ) in symmetric 
nuclear matter (x = 0.5) as function of density u in unit of MeV for (C^* — Cj*)/ 2 = 



for nuclei with x ^ 1/2. This can be done from the results of Friedman et al. by 
noting that our self-energy is nonlinear in x, so 

^(C&,p«po)|*=i/2« 400 &X/&0 MeV (3.41) 

where &o,i are the constants given by Friedman et al. This relation determines the 
coefficient Of* . The result is shown in Table [3] (first three columns) . 

Friedman et al. |38| find the acceptable value to be 61/60 = —0.56 ± 0.82. But 
there is one point which needs to be discussed in interpreting this number in the 
context of our theory. The constant Of* shifts linearly the effective in-medium mass 
of A(1405), with the mass shift being given by 

5m\* = ^2 5E^l(a; = m^* — m B ) (3.42) 

i=a,b 

where 

= - 9 ^\TF K {u)+Tr K {Lo)) 

= -Cl( Pp + Pn ). (3.43) 

For nuclear matter density u — 1 and x = 1/2, the shift is 

5m A * (u, x, y) w [62 - 150.3 x y] MeV (3.44) 

with y = C^f 2 . It seems highly unlikely that the A(1405) will be shifted by hundreds 
of MeV in nuclear matter. This means that y must be of 0(1), and not O(10). For 
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u c 


y = clp 


dAV/dx 


61/60 


F (u) = £ 


F(u) = u 


F(u) = ^ 


50 


125.44Me\/ 


0.314 


2.25 


2.50 


2.97 


40 


64.77Me\/ 


0.162 


2.33 


2.58 


3.08 


30 


4.10MeV 


0.010 


2.42 


2.69 


3.22 


20 


-56.58MeV 


-0.141 


2.54 


2.84 


3.41 


10 


-117.35MeV 


-0.293 


2.71 


3.05 


3.71 


0.41 


-175.43MeV 


-0.439 


2.98 


3.43 


4.28 





-177.92MeV 


-0.445 


2.99 


3.45 


4.32 


-10 


-238.59MeV 


-0.596 


3.60 


4.85 


~ 6.41 



Table 3: Determination of 



from the kaonic atom data]38| and the critical density 
(obtained with the constant so determined) for kaon condensation for various forms 
of symmetry energy F(u) and — CX T )f 2 = 20. y = 0.41 corresponds to no 

A (1405) mass shift in medium at the normal matter density. 



y = 0.41 which corresponds to 61/60 ~ —0.4, there is no shift at normal matter 
density. We believe this is a reasonable value. In fact, y = is also acceptable. It 
would be interesting to measure the shift of A (1405) to fix the constant Cf* more 
precisely although its precise magnitude seems to matter only a little for kaonic atoms 
and as it turns out, negligibly for kaon condensation. 

Let us comment briefly on the role of multi-Fermion Lagrangians. The Weinberg 
counting rule shows that the four- Fermi interactions are suppressed by 0(Q 2 ) relative 
to the terms involving bilinears of Fermi fields. In general n-Fermi interactions will 
be suppressed by the same order relative to (n — 1)-Fermi interactions. In considering 
kaon condensation, what this means in conjunction with the renormalization-group 
flow argument, is that n-Fermi interactions with n > 4 are irrelevant in the RGE 
sense, and hence unimportant for condensation. The situation with the kaonic atom 
data is a bit different. While the strength of the four-Fermi interaction, y, is not 
important (this can be seen in Lee et ol.pl, Table 3), its presence is essential for the 
attraction that seems to be required. This is in contrast to the kaon condensation 
which is driven by the "mass flow" with four-Fermi interactions being irrelevant in 
the RGE sense. 



4 Kaon Condensation 
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4.1 Equation of State and Critical Densities 

We have now all the ingredients needed to calculate the critical density for negatively 
charged kaon condensation in dense nuclear star matter. For this, we will follow the 
procedure given in work of Thorsson, Prakash and Lattimer (TPL)|^T|. As argued by 
Brown, Kubodera and Rho||, we need not consider pions when electrons with high 
chemical potential can trigger condensation through the process e~ — > K~v e . Thus 
we can focus on the spatially uniform condensate 

(K~) = v K e-^. (4.45) 

The energy density e - which is related to the effective potential in the standard way 
- is given by, 

e(u, x, p,v K ) = -Eyuzpo + V(u) + up (l - 2x) 2 S(u) 
5 

- M 2 K - Ilx(/i, U, X)]V 2 K + a n(u, U, X)v n K 

n>2 

+pup x + e e + 9(\p\ - m^e, (4.46) 

where Ef = (pf) /2m B and pf = (3n 2 p /2) 3 are, respectively, Fermi energy and 
momentum at nuclear density. The V(u) is a potential for symmetric nuclear matter 
as described by Prakash et al. [0 which is presumably subsumed in contact four- 
Fermi interactions (and one-pion-exchange - nonlocal - interaction) in the non-strange 
sector as mentioned above. It will affect the equation of state in the condensed phase 
but not the critical density, so we will drop it from now on. The nuclear symmetry 
energy S(u) - also subsumed in four-Fermi interactions in the non-strange sector - 
does play a role as we know from Prakash et al. [f|(| : Protons enter to neutralize the 



charge of condensing K 's making the resulting compact star "nuclear" rather than 
neutron star as one learns in standard astrophysics textbooks. We take the form 



advocated by Prakash et al. [40 



_> 



3 



S(u) = (2i - lj ^^ 0) {u* - F{u)) + S F{u) (4.47) 

where F(u) is the potential contributions to the symmetry energy and Sq ~ 30MeV 
is the bulk symmetry energy parameter. We use three different forms of F(u) ||40|| 



2m 2 

F(u) = u , F(u) = , F(u) = Vu- (4.48) 

1 + u 



The contributions of the filled Fermi seas of electrons and muons are [EI 



l< 



127T 2 

e„ = e„ - up, = ^ ((2t 2 + l)tVWTl - ln(t 2 + - (4-49) 
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clp = 100 


clp = 10 


elf = o 


0.25 


1 


2.25 


3.29 


4.91 


10 


2.25 


3.16 


3.76 


100 


2.18 


2.67 


2.79 



Table 4: Critical density u c in in-medium two-loop chiral perturbation theory for 
F(u) = u. 



where Pf^ = \Jp 2 — is the Fermi momentum and t = pp^/m^. 

The ground-state energy prior to kaon condensation is then obtained by ex- 
tremizing the energy density e with respect to x, p and 



de 
dx 



v K =0 



— 

dp, 







de 



v K =0 



dv 2 



K 







(4.50) 



v K =0 



from which we obtain three equations corresponding, respectively, to beta equilibrium, 
charge neutrality and dispersion relation. The critical density so obtained is given 
for three different F(w)'s in Table |3], and for various ranges of parameters in Table f|. 
The result is 



2 < u c < 4. 



(4.51) 



We note that the largest sensitivity is associated with the part that is not 
controlled by chiral symmetry, namely the density dependence of the symmetry energy 
function F(u). This uncertainty reflects the part of interaction that is not directly 
given by chiral Lagrangians, that is, the part leading to normal nuclear matter. This 
is the major short-coming of our calculation. 

Related to this issue is BR scaling. As we argued, were we able to derive 
nuclear matter from effective chiral Lagrangians, we would have parameters of the 
theory determined at that point reflecting the background around which fluctuations 
are to be made. The BR scaling was proposed in that spirit but with a rather strong 
assumption: That a sigma model governs dynamics in medium as in free space with 
only coupling constants and masses scaled a function of density. Up to date, no 
derivation of this scaling from basic principles has been made. In this sense, we 
might consider it as a conjecture although there is strong support for it from Walecka 
phenomenology in mean field as discussed in ||J]. Suppose we apply BR scaling. The 
only way the procedure can make sense is to apply the scaling argument to the tree 
order terms, but not to the loop corrections. The result of this procedure is significant 
in that the critical density is brought down in an intuitively plausible way to about 
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M& (MeV) 
500 r 




I i i i i i i i i i i i i 

1 2 3 4 5 6 



u = p/p 

Figure 5: Plot of the quantity obtained from the dispersion formula D^ 1 ^, u) = 
vs. the chemical potential /x prior to kaon condensation for g\* = 0.25 and F(u) = u. 
The solid line corresponds to impulse approximation and the dashed lines to the in- 
medium two-loop results for (Cf* — Cj*)/ 2 = 20 and C^f 2 = 20, 10,0 respectively 
from the left. The point at which the chemical potential \x intersects M* K corresponds 
to the critical point. 



u c ~ 2, with very little dependence on parameters, loop corrections and multi- Fermi 
interactions. Thus slightly modified from (|4.51 ), we arrive at the announced result 



2 < u c < 4. (4.52) 



4.2 Irrelevance of A* to Kaon Condensation 

To see which modes are involved in S-wave kaon condensation, we consider the dis- 
persion formula at tree order, 

D- 1 {uj) = uj 2 - M 2 K -Ii(uj). (4.53) 

As shown in Fig. || by the solid line, the kaon "effective mass" is reduced mainly 
by the KN sigma term when there are no A* contributions. If we turn on the A* 
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coupling, there will be additional attractions. However since the effective mass 
lies far from the A*-pole contribution, the resulting magnitude of the attraction is 
small, i.e., the remains nearly unmodified. Furthermore, since the A*-pole is far 
outside of M^, the condensed kaon mode remains the same independently of the A*. 

Summarizing the results, A* may be crucial for understanding the KN scatter- 
ing and kaonic atom data, but is irrelevant to determining the kaon condensation. 
The critical densities for wide ranges of the A* coupling in Table ^ confirm the unim- 
portance of A* contribution. 



5 Kaons on the Hypersphere 



Recently, using the idea of Manton^TJ for simulating density effects, Westerbergp2 
explored S-wave kaon condensation in the bound-state approach to the Skyrme model 
on a 3-sphere. The spatial metric in a hypersphere of radius a is 

ds 2 = a 2 (dp 2 + sin 2 p(d6 2 + sin 2 Odcf) 2 )) (5.54) 

where the possible ranges of three angular coordinates are 

< p, 9<7r, < <p < 2tt. (5.55) 
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Figure 7: Kaon energy (MeV) vs. a [|2 



The baryon number density is given by the inverse volume of the hypersphere 



1 



i F 3 



P 



(5.56) 



2n 2 a 3 2n 2 a 3 

where a = aeF n with F n = 129MeV and e = 5.45. The kaon energy is shown in 
Fig. |7j which also shows that the chiral phase transition occurs at a = a c 
For a > a c , the P-wave and S-wave kaons have different masses, the difference in 
mass representing roughly the mass difference between A(1405) and A(1116). For 
a < a c , the S-wave and P-wave kaons become degenerate, with the kaon condensation 
occurring in the regime 1 < a < 2\[2. 

Solving the equation of state for the electron chemical potential, Westerberg||42|| 
found the critical density to be at a = 1.58 corresponding to p ~ 3.7p Q , Eq. ( |5.56| ). 
This falls within our predicted range (2 ~ 4)po- However an unsatisfactory aspect of 
this result is that the kaon condensation sets in after - and not before as one expects 
- the chiral phase transition. 



6 Discussion 



Introducing A* up to order Q 3 , we obtain the critical density close to that of Kaplan 
and Nelson 0, and confirm that the KN sigma term, as in original approach of 0, is 
essential for kaon condensation. This result is further supported by a renormalization- 
group flow argument as well as by the recent bound-state approach to the Skyrme 
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model. 

Given that kaon condensation occurs at a low enough density as predicted 
here, the Bethe-Brown scenario seems very plausible. However whether or not the 
Bethe-Brown scenario [|J of compact star formation is fully supported by the chiral 
Lagrangian approach will have to await the calculation of the equation of state at 
in-medium two-loop order, which is in progress. Our conjecture is that to the extent 
that our work confirms the original Kaplan-Nelson calculation ||, the compact star 
properties calculated previously at the tree level 0] would come out qualitatively 
unmodified in the higher-order chiral perturbation theory. 
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